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ABSTRACT 

Solitary stars that wander too close to their galactic centres can become tidally disrupted, 
if the tidal forces due to the supermassive black hole (SMBH) residing there overcome the 
self-gravity of the star. If the star is only partially disrupted, so that a fraction survives as a 
self-bound object, this remaining core will experience a net gain in specific orbital energy, 
which translates into a velocity “kick” of up to ^ 10^ km/s. 

In this paper, we present the result of smoothed particle hydrodynamics (SPH) simulations 
of such partial disruptions, and analyse the velocity kick imparted on the surviving core. We 
compare 7=5/3 and 7 = 4/3 polytropes disrupted in both a Newtonian potential, and a 
generalized potential that reproduces most relativistic effects around a Schwarzschild black 
hole either exactly or to excellent precision. For the Newtonian case, we confirm the results 
of previous studies that the kick velocity of the surviving core is virtually independent of 
the ratio of the black hole to stellar mass, and is a function of the impact parameter j3 alone, 
reaching at most the escape velocity of the original star. For a given j3, relativistic effects 
become increasingly important for larger black hole masses. In particular, we find that the kick 
velocity increases with the black hole mass, making larger kicks more common than in the 
Newtonian case, as low-j3 encounters are statistically more likely than high-j3 encounters. 

The analysis of the tidal tensor for the generalized potential shows that our results are 
robust lower limits on the true relativistic kick velocities, and are generally in very good 
agreement with the exact results. 

Key words: black hole physics - galaxies: nuclei - relativistic processes - hydrodynamics - 
methods: numerical. 


1 INTRODUCTION 


The centres of most galaxies are hosts to super massive black holes 
(SMBHs) with masses ranging fr om ^10^ dSecrest et al. 2012| t 
to as much as a few x 10^^ Mq ( van den Bosch et al.||2012 >. The 
SMBH’s mass is comparable to the combined mass of all the stars 
in the nucleus, and together they control the orbital dynamics of 
individual stars. Various mechanisms (the most important being two- 
body scattering; see e.g. | Alexander|2005l l may at times bring one of 
the stars onto a nearly radial orbit that reaches the immediate vicinity 
of the SMBH, with the periapsis dista nce becom ing comparable to 
the tidal radius |Frank 19 ^. On average, this 

happens at a rate of ~ 10“^ yr“^ per galaxy jMagorrian & Tremaine 
|1999[[^ng & Merritt|2004| l. The strength of the encounter, quan¬ 
tified by the impact parameter /3 = r^/r^, will ultimately determine 
how much mass the star loses due to tidal interactions, whether it 
survives (/3 < 1) or is completely ripped apart (/3 > 1) ( |Rees|198^ . 

Recent simulations by [Manukian et al.| l |2013| henceforth 
MGRO; see |Manukian et al.|2014| for errata) showed that in partial 
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disruptions, in which a fraction of the star survives as a self-bound 
object, the remaining core may receive a boost in specific orbital 
energy, corresponding to a velocity “kick” that can reach the surface 
escape velocity of the original star, and may later on be observed 
as a “turbo-velocity star”. The source of this increase in energy has 
been linked to the asymmetry of the two tidal tails created during 
the disruption process, which by conservation of linear momentum 
boosts the velocity of the surviving core. This asymmetry (and the 
velocity kick it induces) appears to be an inherent property of tidal 
disruption events, and has also been observed in tidal disruptions of 
planets ( [Faber et al.|2005|l, white dwarfs ([Cheng & Evans|2013| l and 


neutron stars ( Rosswog et al.|2000|[Kyutoku et al.|2013[ l. 


In this paper we extend the previous results by performing 
smoothed particle hydrodynamics (SPH) simulations of tidal dis¬ 
ruptions of solar-type stars by SMBHs. We use both a completely 
Newtonian approach and one where the orbital dynamics around a 
Schwarzschild black hole is accurately reproduced by a generalized 
potential ( [Tejeda & Rosswog|2Q13[ henceforth TR). Apart from a 
verification of the previous results with different numerical methods, 
our main goal is to quantify to which extent the relativistic effects 
from a Schwarzschild black hole would impact on the final velocit- 
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ies. The study of these events may contribute to the understanding 
of hyper-velocity stars (HVSs; [Hills|1988| l, which were thought to 
result either from the tidal disruption of binary stars, or by scattering 
off the stellar-mass black holes segregated in the Galactic centre 
( [O’Leary & Loeb|2008| ). While these processes can indeed impart 


velocities in excess of 10^ km/s ( [Brown et aL|2005 


Antonini et al. 


|2QTq1 i, most observed HVSs seem to have rather modest velocit¬ 
ies, typically around ^ 400 km/s (e.g. jPalladino et al.|2014[[^eng| 
jet al.|2014||Zhong et al.|2014] l. A number of HVSs have even been 
observed on bound orbits around the Galactic centre, but with suf¬ 
ficiently high velocities so as to constitute a distinct population of 
velocity outliers ( [Brown et al.|2007| l. Tidal disruptions of solitary 
stars may therefore yield sufficiently large kicks to explain many 
of these HVSs, especially if the kicks are significantly enhanced by 
relativistic effects. 


2 METHOD 


In our simulations we use the Newtonian SPH code described in 

detail by[Rosswog et al.[([2009 

1 , with self-gravity computed using 

a binary tree similar to that of 

Benz et al.[([1990[). The tree accur- 


acy parameter (i.e., the opening angle 0 > Hb/Rab which controls 
whether a distant tree node B, of size Hb and located at a distance 
Rab from node A, is allowed to act as a multipole source of gravity 
for node A or needs to be further resolved into its constituents) was 
set to 0 = 0.5, corresponding to an average relative force error of 
<0.1%. All simulations use 10^ SPH particles, unless otherwise 
stated (we verified the results via two runs with 10 ^ particles). 

We model stars as poly tropic fluids with j= 5/3, which 
initially satisfy the Lane-Emden equation for = I Mq and 
r-^ = I Rq. It has long been known that 7=5/3 polytropes are 
disrupted at smaller /3’s than 7=4/3 polytropes, which, being more 
centrally condensed, are able to survive deeper encounters (see e.g. 
[Guillochon & Ramirez-Ruiz[2013[ l. In order to compare our results 
with those of |MGRO who used 7 = 4/3, we also perform a few test 
simulations in which the initial stellar profiles are given by 7 = 4/3 
polytropes, but the fluid, being gas-pressure dominated, reacts to 
dynamical compressions and expansions according to a 7 = 5/3 
equation of state. 

The black hole gravity is modelled with both the Newtonian 
potential (“<Pn’’) and the generalized Newtonian potential (“<Ptr”; 
[Tejeda & Rosswog[2013[ l, 


( 1 ) 

where Es = 2GMbh/c^ is the Schwarzschild radius. The potential 
is based on an expansion of the relativistic equations of motion in 
the low-energy limit, without necessarily implying low-velocities or 
weak field. It has been shown to reproduce essentially all relevant 
orbital properties around a Schwarzschild black hole either exactly 
or to a very high degree of accuracy. 

In the [TR| potential, the specific relativistic orbital energy Etr 
is computed as 


^tr(^,L <P) = - 


GMbh 1 


2 \r-r. 



{r-nY 




r-Vs 


r 


( 2 ) 


The self-bound mass is calculated using the iterative energy-based 
approach described by [Guillochon & Ramirez-RuTz| ( [2013[ Sec. 2 . 2 ), 
with the gravitational self-potential calculated using a fast binary 
tree ( [Gafton & Rosswog[2011[ l. The kick velocity of the self-bound 


object (at infinity) is then computed as 

^kick ~ '\/2(£ — Eq) 1 


( 3 ) 


where Eq is the specific orbital energy of the initial star at the be¬ 
ginning of the simulation, and E is the specific orbital energy of the 
self-bound remnant. This definition does not take into account the 
underlying galactic gravitational potential, and therefore the kick 
velocity is expected to decrease as stars “climb” out of the galactic 
potential. 

Since we are only considering parabolic orbits, Eq is approxim¬ 
ately equal to zero, with small numerical deviations due to the fact 
that only the centre of mass is on a truly parabolic orbit, since we 
impart the same initial orbital velocity to all particles, while starting 
with the star at an initial distance of 5 from the black hole, instead 
of at infinity. 

We simulate encounters with mass ratios q = in the 

range 10 ^ ^^^4x10^, and impact parameters /3 in the range 0.6 ^ 
/3 ^ 0.9 . For the 7 = 4/3 runs we use the same values of /3 and q as 
namely 1.0 to 1.8, and 10^ to 10^, respectively. The values 


MGRO 


for P are chosen so that a self-bound core always survives, while 
q is chosen so that the star is disrupted outside the Schwarzschild 
radius of the SMBH. 

A summary of all the simulations performed for this paper is 
presented in Table 


3 RESULTS 

3.1 Disruption dynamics 

The typical evolution of the stellar fluid during a partial tidal dis¬ 
ruption is shown in Fig. As it approaches periapsis, the star is 
heavily spun-up and distorted, being stretched in the radial direction, 
corresponding to the one positive eigenvalue of the tidal tensor, and 
compressed in the azimuthal and vertical directions (see [Luminet[ 
[& Carter! 1986[ see also Appendix [A| for a discussion on the tidal 
tensor). As the star overfills its Roche lobe it starts to shed mass 
through the Lagrangian points LI and L2, forming a bound and an 
unbound (to the SMBH) tail, respectively. As the star is receding 
from the black hole, the tails and the core stop exchanging energy 
and angular momentum (i.e., the energies become become “frozen- 
in”) and the core recollapses into a self-bound, spherical object. 

Fig. shows the locations of the SPH particles in the time- 
varying Roche potential with respect to the centre of mass of the 
star ('PRoche^ see Appendix [Bjfor a derivation), calculated along the 
radial direction (upper row), and contours of 'PRoche the orbital 
plane (lower row). The formation of the two tidal tails is asymmetric 
from the beginning (the star first overflows its Roche lobe through 
LI), but it is most clearly seen in panel (d 2 ). It is interesting to 
observe that in panel (b 2 ), representative of the star during the actual 
disruption, the Roche potential is not aligned with the star, i.e., the 
points through which the star sheds mass are not always aligned 
with the instantaneous LI and L2. This occurs because the orbital 
time scale of the system is shorter than the dynamical time scale 
on which the fluid can react to the extremely fast-changing Roche 
potential (Torb < %n)- 


3.2 Self-bound mass 

The self-bound mass fraction evolves during the disruption process, 
and for a partial disruption it will decrease from 1 (before disruption, 
when the entire star is self-bound) to the final value mcore/'^^- In 
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Table 1. Overview of the SPH simulations discussed in this paper, grouped into three categories: 59 core 
simulations (1-59) with 7=5/3 and 10^ SPH particles, covering the entire range of q and j 8 discussed in the 
paper; 25 test simulations (60-84) with 7 = 5/3 (of which 24 with 10^ SPH particles, and 1 with 6 = 0.2); 
21 test simulations (85-105) with 7 = 4/3. For each simulation we show the polytropic index 7 , the number 
of SPH particles //part, the ratio q = Mbh/w*, the potential (which is either the Newtonian potential, 
or the generalized Newtonian (TR) potential, <I>tr), and the impact parameter j 8 = n/rp. 


Number 

7 

/^part 



q 

<I> 

e 

i3 

1-7 

5/3 

10 ^ 



lO" 


0.5 

0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90 

8-16 

5/3 

10 ^ 



lO" 

^TR 

0.5 

0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.88, 0.89, 0.90 

17-23 

5/3 

10 ^ 

4 

X 

lO" 


0.5 

0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90 

24-31 

5/3 

10 ^ 

4 

X 

lO" 

^TR 

0.5 

0.60, 0.65, 0.70, 0.75, 0.80, 0.83, 0.84, 0.85 

32-38 

5/3 

10 ^ 



10 ’ 


0.5 

0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90 

39-45 

5/3 

10 ^ 



10 ^ 

^TR 

0.5 

0.60, 0.65, 0.70, 0.75, 0.77, 0.78, 0.80 

46-52 

5/3 

10 ^ 

4 

X 

10 ’ 

<I>N 

0.5 

0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90 

53-59 

5/3 

10 ^ 

4 

X 

10 ^ 

^TR 

0.5 

0.60, 0.65, 0.68, 0.70, 0.71, 0.72, 0.75 

60-66 

5/3 

10 " 



lO" 


0.5 

0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90 

67-73 

5/3 

10 '’ 



lO" 

^TR 

0.5 

0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90 

74-75 

5/3 

10 " 

4 

X 

lO" 

^TR 

0.5 

0.83, 0.84 

76 

5/3 

10 '’ 



10 ^ 


0.5 

0.70 

77 

5/3 

10 " 



10 ’ 

^TR 

0.5 

0.70 

78-83 

5/3 

10 '’ 

4 

X 

10 ’ 

^TR 

0.5 

0.60, 0.65, 0.67, 0.70, 0.71,0.72 

84 

5/3 

10 ^ 



10 ’ 

^TR 

0.2 

0.70 

85 

4/3 

10 '’ 



10 '' 

^TR 

0.5 

1.60 

86 

4/3 

10 " 



lO" 

^TR 

0.5 

1.30 

87-91 

4/3 

10 ^ 



lO" 

^TR 

0.5 

1.45, 1.55, 1.65, 1.75, 1.80 

92-95 

4/3 

10 ^ 



lO" 

<I>N 

0.5 

1.45, 1.55, 1.65, 1.75 

96-100 

4/3 

10 ^ 

4 

X 

lO" 


0.5 

1.10, 1.20, 1.30, 1.70, 1.80 

101-105 

4/3 

10 ^ 

4 

X 

lO" 

^TR 

0.5 

1.10, 1.20, 1.30, 1.65, 1.70 



10^ 


10 ^ 


10-1 ^ 






10 -'' 


10 


-5 


Figure 1. Evolution of the stellar fluid during a typical, parabolic partial tidal disruption. The plot shows a cross-section of the density profile in the orbital 
plane (z = 0). Here, //part = 10^, q = 10^, j8 = 0.7, and = <I>tr. The Schwarzschild radius (black disc), tidal radius (gray disc), geodesic trajectory of the 
centre of mass in Schwarzschild spacetime (dotted red line) and trajectory of the centre of mass obtained with the TR potential (dashed black line) are shown to 
scale. Due to the spatial scales involved, we show the stellar debris magnified by a factor of 50 before periapsis passage {t ^ 0) and magnified by a factor of 10 
afterwards (f > 0). Note that due to relativistic periapsis shift the orbit is not a parabola. 
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Figure 2. A typical, time-varying Roche potential in a partial tidal disruption exhibits a number of stages, shown here as snapshots at: (a) t — 3.5 h before 
periapsis passage, i.e. at the beginning of the simulation, when the star is self-bound; (b) t min, just as the first particles exit the Roche lobe of the star and 
start forming the tidal tails (bound tail first); (c)t ^4 h, as the energy distribution of the debris starts to freeze and the core and tails approach their final masses; 
(d) t ^4 days after the disruption. This simulation used Apart = 10^, ^ = 10^, j8 = 0.65, and = <I>n. Upper row. 'RRoche (^)5 where r is the distance to the 
SMBH. The values of 'RRoche are given in code units and - to simplify axis labelling - are offset by the values shown in the upper left comer of the panels. 
Lower row. Contours of 'RRoche (-^7 j) in the orbital plane, with the particles overplotted and coloured according to the value of the potential. All coordinate axes 
use a global Cartesian coordinate system, with the SMBH always located at (0,0,0), and the stellar fluid moving along a parabolic orbit around it. In the lower 
row we are essentially “zooming in” on the SPH particle distribution as it first approaches, and then recedes from the SMBH. 




Figure 3. Left panel. Evolution of the self-bound mass fraction in simulations using Newtonian (dashed lines) and |TR| (solid lines) potentials. Right panel. Final 
self-bound mass fraction in Newtonian (dashed lines) and |TR| (solid lines) simulations. The points corresponding to Newtonian runs are esentially overlapping 
for all values of j8, while those from pseudo-relativistic simulations show a strong dependence on q: the larger q is, the smaller the critical j8 necessary for 
complete disruption is. 


Fig. (left panel) we present the time-evolution of mcore to the 
^ = 10^ simulations. The disruption is stronger in relativistic en¬ 
counters, with the deviation of the self-bound mass fraction from the 
Newtonian case increasing with /3, from ^ few percent (/3 = 0.6) to 
^100 percent (/3 =0.9). 

Fig.[^(right panel) shows the final self-bound mass fraction as a 
function of /3, for various ratios q and for both potentials. We observe 


that with the Newtonian potential stars are partially disrupted in 
the range 0.6 < /3 < 0.9 regardless of q. This result is general for 
all Newtonian disruptions of 7 = 5/3 polytropes, and agrees with 
the numerical findings of e.g. |Guillochon & Ramirez-Ruiz| ( [T013| l 
(/3d necessary for complete disruption equal to 0.9). On the other 


hand, for any given /3 the discrepancy between Newtonian and TR 


simulations increases drastically for larger black hole masses, as r^ 
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becomes comparable to and relativistic terms in the tidal tensor 
close to periapsis can no longer be ignored (see Appendix [A]). 


3.3 Kick velocity 

The kick velocity is computed from the increase in specific orbital 
energy (e — Eq) of the self-bound core using Eq. In Fig. we 
present v^ick both as a function of time for the ^ = 10 ^ simulations 
(left panel), and its final value at the end of the simulation (right 
panel). The kick velocities are always scaled to the surface escape 
velocity of the original star, sJIGMq/Rq '^611 km/s, as they are 
expected to be comparable on theoretical grounds (see|Manukian| 
let al.|2013| ). 

We initially place the star at 5 and observe that numerical de¬ 
viations from Co = 0 (i-e., parabolic trajectory) of the initial specific 
energy increases with q, and is larger for the |TR| potential. For the 
most extreme case {q = Ax 10 ^, /3 = 0 . 9 , the initial energy 
was of the order of 10 “^ (which is to be compared with the 
typical “kick” specific energy, ~ 10 “^ c^), while for the other simu¬ 
lations it was between one and three orders of magnitude smaller 
(for the |TR| and Newtonian potentials, respectively). Nevertheless, 
there is a clear distinction between the Newtonian runs, where the 
kick is fairly independent of q, and the |TR| runs, where the curves 
exhibit an asymptotic behaviour limited by progressively smaller 
/3’s with larger ^’s. 

The time evolution of e and - consequently - of v^ick exhibits 
significant oscillations during the actual disruption (~ hours from 
periapsis passage) that appear as wiggles in Fig.|^ The explanation 
is likely related to the complicated hydrodynamic effects that take 
place in that short period of time. Due to the extreme compression 
of the star, shock waves travel throughout the star and transfer 
significant energy and angular momentum between the particles, 
and until the energies become frozen-in, the average specific energy 
of the bound core keeps oscillating. 

Fig. 1^ shows the kick velocity as a function of the self-bound 
mass fraction mcore- We notice that in contrast to the findings of 
|MGRO| in some simulations the kick velocity exceeds the escape 
velocity of the initial star, but only when the value of mcore is suf¬ 
ficiently small. Indeed, Vkick('Wcore) seems to be a monotonic func¬ 
tion that asymptotically approaches 0 for mcore ^ and +00 for 
^core 0- 

Fig. shows the kick velocity as a function of the mass dif¬ 
ference between the two tidal tails. Am 12 (Fig. left and right 
panels, for Newtonian and |TR| potentials, respectively). The two 
plots exhibit similar behaviours (v^ick increases with Ami 2 ), but in 
the relativistic simulations the degeneracy in q is broken (i.e., the 
data points for a given /3 are not clustered together irrespective of 
q), since the relativistic kicks are sensitive to both /3 and q. 

We also show v^ick as a function of /3 for the 7=4/3 simu¬ 
lations (Fig.[^, together with the fit line given by |MGRO| Th ese 
simulations show reasonable similarity to the results ^MGRO| and 
- as expected - for such small values of q ( 10 ^ to 10 ^) there is little 
difference between Newtonian and |TR| simulations. The data points 
from our simulations give slightly smaller kicks than the fit, but the 
general trend (and critical /3) are nevertheless recovered. 


3.4 Error estimation and resolution dependence 

Since in this paper we are concerned with parabolic encounters, 
we have calculated the |TR| acceleration for particles on parabolic 
orbits. It turns out that the TR potential reproduces parabolic orbits 



mcore/m,, 


Figure 5. The kick velocity Vkick of the self-bound remnant shows a remark¬ 
able dependency on the the self-bound mass fraction in both Newtonian 
(blue points) and |TR| (red points) simulations. We fitted a truncated power 
law of the form y = (1 — , with the fit parameters A = 0.634726, 

a = 0.196598, j8 = 0.882387. The shaded gray area around the black fit line 
represents the 1 cr deviation from the fit. Darker points correspond to higher 
values of j8. 


0.8 - 




0.2 - 





f 





• ■ 
u' 

.A. 

1 




i • 


_. 1 

-A'' 

ir.:.;. . 




,. - I 

i 

_1 

1 w 

■ 

i_i 

1_1 

i_1 

1_ 


1 1.2 1.4 1.6 

Impact parameter j8 

Figure 7. The results of our SPH simulatio ns of 7 = 4/3 polytropes, using 
^^part = 10^ in Newtonian (black circles) and 1 


TR 


(red squares) potentials, and 
using Apart = 10^ with the TR potential (red triangles), as compared to the 
fit line given by MGRO for their 7 = 4/3 simulations. Our runs use the same 
values for q as in their paper (10^ to 10^), and for this reason relativistic 
effects are small. 


exactly (see |Tejeda & Rosswog||2013[ in our simulations, these 
would be the orbits of the centre of mass). In Fig. we show the 
tidal deformation experienced by test particles placed on a ring 
of radius ARq with the centre following a parabolic orbit with 
/3 = 1.6 around a 10^ Mq black hole. The deformation of the 
ring provides a visual representation of the effects of the tidal 
tensor, and by eye it appears to be identical in Schwarzschild 
and In Appendix [A| we calculate in detail the components 
and the eigenvalues of the tidal tensor in Schwarzschild and both 
Newtonian and |TR| potentials, and show that the eigenvalues A 2 
and A 3 (corresponding to compression in the orbital plane and 
perpendicular to it, respectively) provided by the |TR| potential 
for parabolic orbits are the same as in Schwarzschild (with A 3 
having an identical expression, and A 2 having a slightly different 
expression that still gives the same result within machine precision 
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Figure 4. Left panel. Evolution of the kick velocity of the self-bound remnant in Newtonian (dashed lines) and |TR| (solid lines) simulations. The kick velocity is 
normalized by the surface escape velocity of a solar-type star, 617 km/s. The lines show a moving average of the data points in order to smooth out fluctuations 
during the disruption process (before the surviving remnant becom es relaxed, i.e. approximately during ^ = 0 and f = 1). Right panel. Kick velocity of the 
self-bound remnant at infinity, in Newtonian (dashed lines) and |TR| (solid lines) simulations. The solid black line {q = Ax 10^, = ^tr) stands out, as for low 

values of j8 it suprisingly falls both below the runs with smaller q, and below the Newtonian run. We have repeated this set of simulations with higher resolution 
(10^ SPH particles), but the results were very similar. Since these are highly relativistic encounters, further studies with even higher resolution and an exact 
relativistic treatment of the black hole gravity are probably required in order to arrive at a definitive answer.. 
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Figure 6. Kick velocity of the self-bound remnant as a function of the mass difference Amn between the tidal tails in Newtonian {left panel) and |TR| (ng/?r 
panel) simulations. The value of j8 is colour-coded. The relativistic results do not exhibit such a strong relation between j8, Ami 2 , and Vkick as in the case of the 
Newtonian potential, but they still obey the general trend (vkick increases with Ami 2 ) and even the shape of the function. 


for our range of parameters), while the maximum relative error in 
X\ (corresponding to expansion in the radial direction) is 6.4%, 
occurring at periapsis for a very narrow range of /3’s. In general, 
the average error in X\ at periapsis is of order ~ a percent, and 
quickly drops farther away from the black hole. We also show that 
both the Newtonian and the |TR| potentials always underestimate 
the values of X\ and A 2 , which means that all disruptions will be 
slightly stronger in the Schwarzschild spacetime, leaving a smaller 
self-bound core and resulting in even larger velocity kicks. There¬ 
fore, our results are robust lower limits on the true relativistic effects. 

We have also performed a number of test simulations with higher 
resolution. For ^ = 10^, we ran all simulations with both 10^ and 
10^ SPH particles; a comparison of the resulting kick velocities is 
presented in Fig.[^ The most striking observation is that resolution 
is crucial at low /3 (~ 0.6), where the mass asymmetry between the 
two tidal tails, which drives the kick, is extremely small, < 10“^ m^. 


This mass loss must be properly resolved in order for the kick effect 
to be accurately captured; we observe the low-resolution simulations 
greatly underestimate the kick for /3 = 0.6, due to more energy 
being deposited into the oscillation modes of the star. Resolution 
is also important, though not as crucial, at high /3 (~ 0.9), where 
the surviving core itself is of the order of ^ 10“^ m^, and must also 
be properly resolved; we observe that here the discrepancy between 
On and Otr is less extreme for the higher resolution simulations, 
but this is simply because Otr leaves a small surviving core, which 
is not resolved by the low-resolution simulations. For intermediate 
values of /3, where both the core and the mass difference of the tails 
are a generous fraction of the initial stellar mass, lower-resolution 
simulations agree to a reasonable accuracy ('^ few percent) with the 
results obtained with the 10^ particles simulations. 

For q = 10^, we ran the simulation with /3 = 0.7 with vari¬ 
ous numbers of SPH particles A^part and tree opening angles 0: for 

^part = 10^ and A^part = 10^ (both with 0 = 0.5), and for Otr, 
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Figure 8. Visual representation of the tidal deformation experienced by a ring of particles of radius ARq (representing a star disrupted on a parabolic orbit with 
j8 = 1.6 by a SMBH with q = 10^) in Schwarzschild spacetime (red), the |TR| potential (blue), and the Newtonian potential (green). The rings in Schwarzschild 
spacetime and the TR potential are indistinguishable from each other by eye, which is an indication that <I>tr reproduces the Schwarzschild tidal tensor to a very 
good accuracy (see Appendix [A| for a further discussion). Unlike in the Newtonian case, where the shape of the rings would not change significantly after 
they exit the tidal radius, the other two rings become increasingly deformed due to the different periapsis shifts of their component particles (the TR potential 
reproduces the periapsis shift exactly). 




Figure 9. Left panel. Kick velocity as a function of impact parameter j8, obtained with lower (Apart = 10^; black lines) and higher (Apart = 10^; red lines) 
resolution simulations, representing both <I>n (dashed lines) and <I>tr (solid lines) encounters. Resolution appears to be critical at low j8 (~ 0.6), where the 
mass asymetry of the two tails must be properly resolved in order to extract an accurate kick velocity (nevertheless, at both resolutions there is a clear trend of 
higher kick velocities in the relativistic case, compared to the Newtonian case). For the central part of the range, simulations with the two resolutions agree to 
within ~ 1 percent. At large j8, close to the critical impact parameter (~ 0.9), the two results agree for but for the TR encounter the 10^ particles 

simulation results in no surviving core, while the 10^ particles simulation still yields a small self-bound mass that acquires all the energy, resulting in a very large 
kick. Right panel. Time evolution of the kick velocity fox 10^, = <Ftr, and Apart = 10^ (dashed lines) and Apart = 10^ (solid lines). The different colours 

represent different values of j8. The best agreement between the kick velocities is achieved in the region of intermediate j8, where both the mass asymmetry of 
the two tails and the surviving core are resolved in the lower-resolution runs. Since a significant fraction of the star survives for low j8’s, it is very expensive to 
run these simulations with high resolution, and therefore they were stopped at an earlier physical time than the high j8 simulations. 
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^part = 10^ = 0.2, 0.5), and A^part = 10® (0 = 0.5). We observe 

that improving the force accuracy above 0 = 0.5 does not have a 
significant impact on the final kick velocity (at most ~ one percent), 
while increasing the number of particles ten times does change it 
by up to ^ 10 percent, for both potentials. In addition, the evolu¬ 
tion of Ekick during the disruption process is highly dependent on 
the resolution, with larger wiggles in the high-resolution simula¬ 
tions, presumably due to the fact that shocks are better resolved 
and therefore less dissipative, and there is significantly more energy 
transfer between the particles. Still, the difference between the kick 
velocities of the Newtonian and |TR| potentials is consistently and 
significantly larger than the variations that appear when changing 
the accuracy parameters for the same potential. 

For 7 = 4/3, we ran two additional simulations with A^part = 
10^ and 0 = <Ftr, using q = 10"^, /3 = 1.6 and q = 10^, /3 = 1.3. 
The kicks obtained in these higher-resolution simulations are shown 
in Fig. 1^ with red squares, and are in good agreement with both the 
fit given by MGRO and the rest of the points from our simulations. 


4 DISCUSSION 


We have confirmed by a different numerical method (SPH vs. the 
adaptive mesh refinement code FLASH that is primarily used by 
|MGRO| ) that all Newtonian simulations with the same impact para¬ 
meter P produce similar self-bound remnants (in mass, specific 
energy, and kick velocity), confirming earlier conclusions by |Ma-| 
|nukian et al.| ( |2Q13} that Vkick is virtually independent of the mass 
ratio q. 

We compared these calculations with simulations that use a 
generalized potential that accurately captures the dynamics around 
a Schwarzschild black hole. Unlike Newtonian tidal disruptions, we 
observe that relativistic disruptions are no longer fully described 
by the parameter P . Instead, relativistic effects related to the tidal 
tensor become important when the periapsis distance is comparable 
to the Schwarzschild radius, i.e. they depend on the ratio 


0E 






1 ' m* 

( ^ "l 

s' 

KReJ 

' Me 



2/3 


• (4) 


For a given mass ratio q, this can be interpreted as a depend¬ 
ence on P (e.g., for ^ = 10^, relativistic effects become important 
when P >5, as noticed by e.g. [Laguna et al.||1993| l, which will 
however change with q (for q = 10^ relativistic effects are extremely 
important even for /3 < 1), since the gravitational radius and the 
tidal radius have different dependencies on M^h. In general, we ob¬ 
serve that relativistic effects can be ignored for 0 < 10“^ (where 
the relativistic and the Newtonian tidal tensor are essentially the 
same, and relativistic effects such as periapsis shift are negligible on 
the time scale of the disruption), but tend to dominate the outcome 
of the encounter for 0 > 0.1. These thresholds can also be observed 
in Fig. |A1| in the Appendix, where we see that in a typical tidal 
disruption the eigenvalues of the tidal tensor are virtually identical 
in Newton and Schwarzschild for r > 30 rg, and start to differ signi¬ 
ficantly from each other for r < lOrg (the exact values will of course 
depend on the parameters of the encounter). 

We therefore expect relativistic effects on the kick velocity of 
turbo-velocity stars to be significant for larger black hole masses. In 
particular, since these objects will only result from partial disrup¬ 
tions (0.6 < /3 < 0.9, with the more pronounced velocity kicks at 
the higher end of this range), relativistic effects should dominate for 


0 > 0.1, or, from Eq. (|^, for 


^>4.2x10^ 


Me. 




\Re 


3/2 


(5) 


This means that for a solar-type star disrupted in our Galactic centre 
relativistic effe cts m ay be important. In Fig. we observe that for 


^ = 4 X 10” the TR potential will yield the same kick velocity as the 


Newtonian potential at a ^ 5% smaller p. Since tidal disruption rates 
scale with ~ /3“^ we estimate that approximately 10% more stars 
will have a given kick velocity compared to a Newtonian estimation. 

Previous estimations by |Kesden| ( |20 1 2a|b'] t predict that the spin 
of the black hole may alter the spread in energy by up to a factor of 
^ 2. If one were to anticipate the relativistic effects due to the black 
hole spin, one would therefore expect a maximum of ^ 40% increase 
in the kick velocity, depending on the spin and orbit orientation, 
though a methodical study of such effects is left for subsequent 
investigations. 

To conclude, we have found that the critical P necessary for 
the disruption of the star is highly dependent on q as long as 0 > 
10“^, with the star being completely disrupted at a /3 of around 0.9 
(q = 10^), 0.85 (^ = 4 X 10^), 0.8 (q = 10^), 0.72 (^ = 4 x 10^), 
due to the proximity to the event horizon (0 ^ 0.04, 0.09, 0.16, and 
0.36, respectively). This implies that the higher q is, the smaller the 
span of P's in which partial disruption will occur, and the steeper 
the dependency of the surviving core’s mass, mcore on P is (Fig.[^. 
Since we have also shown that there is a very clear, monotonic 
relation between mcore and the kick velocity Vkick imparted to the 
core (Fig. [^, we conclude that - in Schwarzschild spacetime - 
heavier black holes are able to impart larger kick velocities without 
requiring very deep encounters. 
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APPENDIX A: DERIVATION OF THE TIDAL 
GRAVITATIONAL FIELD FOR THE TR POTENTIAL 


In this appendix we provide explicit expressions for the tidal tensor 
corresponding to the |TR| potential and quantify its departure from the 
exact relativistic result in Schwarzschild spacetime. The acceleration 
exerted on a given test particle under the TR potential is given by 
( [Tejeda & Rosswog|2013| ) 


GMbh^,- ^ rsXjr 3 r^Xj (p^ 

V rJ r(r-rs) 2 r 


(Al) 


where Xf = {v, y, z} and r, r and (p should be taken as implicit 
functions of the Cartesian coordinates satisfying 

(A2) 

rr = xx-\-yy-\-zz, (A3) 

= {xy-yxf + {xz-zxf + {zy-yzf. (A4) 


We now consider a star approaching the central black hole. The 
centre of mass of the star is located at xq and follows, to a very 
good approximation, the trajectory of a free-falling test particle. 
We can compute the tidal forces acting on a fluid element located 


at a generic position x within the star by taking the Taylor series 
expansion of Eq. (AT) around xo, i.e. 


Xi=fi{x,x)z^ + ^ 




, . s, dfi 


(j^oAo) 


(A5) 


where the Einstein summation convention is used. Next we substitute 
^ =x — Xq into Eq. ( |A5j ) and find that, to first order in ^ and ^, the 
acceleration acting on a fluid element due to the central black hole 
as seen from the centre of mass of the star is given by 


Q 






(ji^oAo) 


(A 6 ) 


with the tidal tensors Qj and Qj given by 


dfi 

dXj 


GMbh 



iS — 

f3r-5n 

\ Xixf 
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r / 
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\ r-r. 

) . 

rsXiXj 

Sij 

_ 1 

Or- 

IrA rxj' 


d(r-n) 



r, J rxj 



(A7) 


dfi 

dXj 



where 5ij is the Kronecker delta symbol and + 3 )^ +z^. 

Now we would like to compare the results in Eqs. ( |A7j l and 
( |A 8 | l with the corresponding exact relativistic expressions. It is cus¬ 
tomary to express the tidal tensor in Schwarzschild spacetime by 
adopting the so-called Fermi normal coordinates centred on the 
approaching star (see e.g. |Marck|1983[|Brassart & Luminet|2010| l. 
Nevertheless, for consistency with the approach that we have adop¬ 
ted here, we need to calculate the tidal field using Schwarzschild 
coordinates (i.e. in the global reference frame centred on the black 
hole). The starting point is then the full-relativistic expression for 
the acceleration acting on a test particle in Schwarzschild spacetime 


.. ^ GMbhX/ nxj r r^Xir"^ _ r^Xj 

P \ r) r{r-n) 2{r-r,)d r 

(A9) 

Following the same steps leading to Eq. ( |A 6 j ), it is found that the 
tidal tensors are now given by 



+ 


GMbh 


• - 1 

(3r-4rA 

XiXj- 

r3 V 

r J 
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K r-rs J 

1 ^2 

r,XiXj 

Sij + 

rxf 

/3r-2rZ 

\ 

d{r-r^) 

rxf 

V , 

/ rxj 



J r\r-r,) ^4 v y 


(AlO) 


(All) 


where the superscript ‘S’ is used to indicate that the quantity has 
been calculated in Schwarzschild spacetime. 

The eigenvalues of the tidal tensor Qj give pertinent physical 
information about the amount of compression or expansion that the 
stellar matter experiences due to the black hole along the direction 
of the principal axes of the tidal tensor. These eigenvalues can be 
found by following the standard procedure of diagonalizing Qj. In 
particular, if we choose a global reference frame XYZ such that 
the trajectory followed by the centre of mass of the star coincides 
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Figure Al. Comparison of the eigenvalues of the tidal tensor as obtained in 
Schwarzschild spacetime (solid lines), for the |TR| potential (dashed lines), 
and for the Newtonian potential (dotted lines). In this case we have taken 
a parabolic encounter with h = Sr^c, corresponding to a j8 ~ 1 encounter 
between a solar-type star with a lO^M© central black hole. 


with the XY plane (which is always possible due to the spherical 
symmetry of the present problem), the corresponding eigenvalues 
are given by 


^ 1,2 = 


Q 


1 + C22 ± \/(Ql +<^22)^ - 4 (CiiC 22 -C12C12) 


h =C33. 


(A12) 


Equivalent expressions are found for 2.^ by substituting Cf- instead 
of Cij into Eq. ( |A12| ). 

In Eig. |Al| we compare the eigenvalues in Eq. ( |A12| ) with the 
corresponding relativistic values Xf for a parabolic trajectory with 
h = Sr^c (which represents a /3 ~ 1 encounter between a solar-type 
star and a 10^ Mq black hole). As a reference to illustrate the import¬ 
ance of relativistic effects, we have also plotted the corresponding 
Newtonian values ( [Brassart & Luminet|2008| l 


^ ^ _ 2 GMbh 

^1 “ 


a A _ a A _ ^^bh 
A2 - A3 - ^ 


(A13) 


In this figure and for any other parabolic encounter, the eigenvalue 
A 3 coincides exactly with the corresponding relativistic result. More 
specifically, for a parabolic trajectory with angular momentum h, it 
is found that 


A3 = A| 



(A14) 


This means that tidal compression along the vertical direction is 
reproduced exactly by the |TR| potential for parabolic encounters. On 
the other hand, it is also apparent from Eig. |Al| that the eigenvalues 
Al and A 2 provide a good approximation to the exact relativistic 
values. 

In a practical simulation, the maximum departure of the tidal 
field corresponding to the |TR| potential from the relativistic result 
is reached at the periapsis of the incoming trajectory. In Eig. |A2| 
we show the contours of the maximum relative error for Ai (i.e. 
(Aj^ — Al) / Aj^) as a function of the impact parameter /3 and the mass 
ratio q = From this figure we see that expansion due to the 

tidal tensor is reproduced by the |TR| potential with a precision better 
than 6.4%. Finally, the relative error for A 2 at periapsis was found 
to be consistently zero within machine precision (p:^ 10 “^^) for the 
parameter values shown in Fig. |A2[ This appears to contradict the 


Figure A2. Contours of the maximum relative error with which the eigen¬ 
value of the tidal tensor Ai (Eq. EH reproduces the exact relativistic value 
Af as a function of the impact parameter j8 and the mass ratio q = Mbh/m*. 


fact that Aj^ deviates from A^ at r ~ 2r^ in Fig. 


Al 


; at pericentre, which in Fig. |Al| is at 


maximum error occurs « 
and that the curves shown there will change with /3. 


but note that the 
^ 24n, 


APPENDIX B: THE ROCHE POTENTIAL EOR 
NON-CIRCULAR ORBITS 


The Roche lobe is defined in the context of the so-called restricted 
three-body problem, where a test particle moves in the potential of 
two orbiting masses. It is bounded by the critical equipotential sur¬ 
face of the effective potential d>eff(r), incorporating inertial forces 
in the coordinate frame that is comoving and corotating with the star 
( |Sepinsky et al.|20()7| l. In such a non-inertial frame, the expression 
for the effective potential reads: 

1 . 

<l>eff ('■)=<!>* (r)+ <I>bh V* <l>bh (/■*)- z I ® ('■) X (r - r*) I . 

(Bl) 

Here, r* is the current instantaneous position of the centre of mass 
(CoM) of the star, V* is the gradient with respect to r*, (o{r) is 
the angular velocity at the position r. The first two terms in this 
expression, and d>bh(^), represent the gravitational potentials 
of the star and the black hole, respectively; in our simulations, the 
former is computed with the tree, while the latter can be either 
d>N or Otr. The third term appears because our reference frame 
is comoving with the star. It produces uniform acceleration, equal 
and opposite to the one of the CoM of the star. The last term is the 
centrifugal potential due to stellar rotation. 

If we use point particle Newtonian potentials for the black hole 
and the star, the effective potential becomes: 


Gm* GMbh GMbh 

'^'effW =- 

r r ri 


{r-n)- -\(0{r)yi{r-r^)\^. 


(B2) 

The choice of angular velocity ft) in the effective potential de¬ 
pends on how the star is rotating and is not very clearly defined 
in the post-disruption phase. Therefore, in Fig.[^ for the contours 
of d>eff(r) (dashed black lines) we adopt an average angular velo¬ 
city (ft)) = (where a is the particle index), while for the 

colours of individual particles we use the values of the individual 
angular velocities ft)^. 
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